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Dilepton invariant-mass spectra for heavy-ion collisions at SIS 18 and BEVALAC energies are 
calculated using a coarse-grained time evolution from the Ultra-relativistic Quantum Molecular Dy¬ 
namics (UrQMD) model. The coarse-graining of the microscopic simulations enables to calculate 
thermal dilepton emission rates by application of in-medium spectral functions from equilibrium 
quantum-field theoretical calculations. The results show that extremely high baryon chemical po¬ 
tentials dominate the evolution of the created hot and dense fireball. Consequently, a significant 
modification of the p spectral shape becomes visible in the dilepton invariant-mass spectrum, re¬ 
sulting in an enhancement in the low-mass region M ee = 200 to 600MeV/c 2 . This enhancement, 
mainly caused by baryonic effects on the p spectral shape, can fully describe the experimentally ob¬ 
served excess above the hadronic cocktail contributions in Ar+KCl (Ui a b = 1.76 AGeV) reactions as 
measured by the HADES collaboration and also gives a good explanation of the older DLS Ca+Ca 
(_E lab = 1.04 AGeV) data. For the larger Au+Au (i?iab = 1.23 AGeV) system, we predict an even 
stronger excess from our calculations. A systematic comparison of the results for different system 
sizes from C+C to Au+Au shows that the thermal dilepton yield increases stronger (oc A 4 / 3 ) than 
the hadronic background contributions, which scale with A, due to its sensitivity on the time evolu¬ 
tion of the reaction. We stress that the findings of the present work are consistent with our previous 
coarse-graining results for the NA60 measurements at top SPS energy. We argue that it is possible 
to describe the dilepton results from SIS 18 up to SPS energies by considering the modifications of 
the p spectral function inside a hot and dense medium within the same model. 
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I. INTRODUCTION 

The study of dilepton production has for long been pro¬ 
posed as a good method to probe the change of hadronic 
properties in the hot and dense matter created in heavy- 
ion collisions and also as a possible observable for the 
creation of a deconfined phase at sufficiently high colli¬ 
sion energies UHH. In contrast to the vacuum situation, 
a hadron can not only decay in a hot and dense medium 
but also interact with the constituents of the medium in 
scattering processes and resonance excitation. Many the¬ 
oretical efforts have been undertaken over the last years 
to gain a better understanding of the in-medium proper¬ 
ties of vector mesons EHB]. The behavior of hadrons is of 
immanent interest for a full understanding of the phase 
structure given by Quantum Chromodynamics (QCD). 
One much-discussed aspect here is the change of the sym¬ 
metries of QCD which is expected when going from the 
vacuum to finite temperature and finite baryochemical 
potential, especially the predicted restoration of chiral 
symmetry mm- Unlike all hadronic observables, which 
only provide information on the final freeze-out stage of 
the system, dileptons are not subject to strong inter¬ 
actions and consequently escape the fireball unscathed. 
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However, this also means that lepton pairs from all stages 
of the reaction will reach the detector. Especially for a 
theoretical description this is a big challenge, as it de¬ 
mands a realistic description of the whole space-time evo¬ 
lution of the heavy-ion reaction and taking the various 
dilepton sources into account. 

On the experimental side several groups have under¬ 
taken the challenging task to measure dilepton spectra in 
heavy-ion collisions and thereby constrain the theoreti¬ 
cal predictions. At SPS energies the NA60 Collaboration 
was able to measure the p in-mediunr spectral function 
for the first time, thanks to the high precision of the 
measurement m- The results were in line with previous 
CERES results |l2j and found an excess in the invariant 
mass range from 0.2 to 0.6GeV/c 2 . This excess can be 
explained by a broadening of the p meson inside the hot 
and dense medium with small mass shifts p~3l - fT6] . At 
RHIC, these investigations were extended to even higher 
collision energies with basically the same results except 
for less dominant baryonic effects and a larger fraction of 
dileptons stemming from the QGP puns]. 

Still more challenging is the interpretation of the dilep¬ 
ton measurements, which were performed in the low- 
energy regime at SIS 18 and BEVALAC. For collision 
energies around E\ a ^ ~ 1 — 2AGeV, which will be in 
the focus of the present study, the DLS Collaboration 
measured a large excess beyond the results of theoreti¬ 
cal microscopic calculations several years ago [Tj3]. This 
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disagreement between experiment and theory was called 
the “DLS puzzle”. In consequence, it triggered further 
experimental and theoretical investigations. Recently, 
the HADES experiment confirmed the former DLS re¬ 
sults with a higher precision | EUH25] , Although the the¬ 
oretical microscopic models have been largely improved 
and extended since this time [MH23] . a full and unam¬ 
biguous description of the data has not yet been found. 
A satisfying answer is complicated by the fact that at 
such low energies a large number of processes contributes 
to the dilepton production, for which many parameters 
(like cross-sections, branching ratios, etc.) are not well- 
known. In addition interference effects are posing serious 
problems for transport Monte-Carlo simulations based 
on the evolution of phase-space densities. Here future 
measurements, for example in pion induced reactions as 
conducted by HADES, could give better constraints for 
the various parameters and reduce the uncertainty of the 
different contributions [301. In any case, a full descrip¬ 
tion of in-medium effects via off-shell dynamics or multi¬ 
particle interactions at high densities remains a difficult 
task, although some investigations on these issues have 
been conducted successfully [25JI3TH37]. 

Besides the microscopic transport models, there ex¬ 
ist also macroscopic approaches describing the evolution 
of the heavy-ion collision in terms of its thermodynamic 
properties. At high collision energies, i.e. at SPS or 
RHIC, a thermal calculation of dilepton emission is often 
applied, where a fireball expansion or a hydrodynamic 
calculation is used to model the bulk evolution of the 
system, while the dilepton emission is calculated using 
spectral functions at a given T and pb H3MI51I3S0DI. 
But this approach works only if the collision energy is 
high enough. The application at ss 1-2 AGeV is 
hardly reliable. 

However, there are no reasons why a macroscopic de¬ 
scription of the reaction dynamics should not be possible 
at SIS 18 and BEVALAC energies, provided one can ex¬ 
tract realistic values of energy and baryon density and, in 
consequence, temperature and baryocliemical potential. 
On the contrary, due the the expected high values of the 
baryon chemical potential, a study of the thermal proper¬ 
ties of the system and the influence on the spectral shape 
of vector mesons might be very instructive. In the present 
work we argue that it is possible to obtain realistic values 
of T and pb from microscopic calculations, provided one 
uses a large ensemble of events and averages over them 
(i.e. one ”coarse-grains” the results) to obtain a locally 
smooth phase-space distribution. This ansatz, which was 
first presented and applied for the calculation of dilep¬ 
ton and photon spectra in SB. constitutes a compromise 
between (non-equilibrium) microscopic transport simula¬ 
tions and the calculation of clilepton emission with (equi¬ 
librium) spectral functions. 

The same model was described in detail in 16 and 
already successfully applied to investigate dilepton pro¬ 
duction at SPS energies [HIS3]. The ansatz is used ba¬ 
sically unchanged for the present low-energy study. The 


only important extension is the implementation of the 
in-medium spectral function for the to meson. While for 
the NA60 results the cocktail contribution of the u was 
already subtracted from the thermal dilepton spectra, for 
comparing our calculations to the experimental HADES 
and DLS results a full description of the to contribution 
is required. Furthermore, the high baryon densities and 
slow evolution of the system increase the significance of 
its medium modifications at the energies considered here. 

Note that we focus the present investigation on the 
larger systems Ar+KCl and Au+Au, where local ther- 
malization can be assumed due to the size of the hot 
and dense fireball. For the smaller C+C system, which 
was also studied by HADES and DLS (in nearly mini¬ 
mum bias reactions), one finds that the average number 
of NN-collisions is so small that the assumption of a local 
thermalization is questionable. Furthermore, the exper¬ 
imental results indicate that the C+C dilepton spectra 
can be interpreted as the superposition of the underly¬ 
ing p+p and p+n collisions without any significant in¬ 
medium effects [HI- (However, when studying the effect 
of the system size on the dilepton production we will 
also consider central C+C collisions later on for com¬ 
pleteness.) 

This paper is organized as follows. Section |TI| gives an 
introduction to the model and an overview of the different 
contributions considered for the calculation. Section 1III1 
then presents the results for the thermodynamic evolu¬ 
tion of the reaction (Sec. 


Ill A) and the resulting dilepton 


spectra (Sec. Ill B), furthermore the effect of the system 
size and fireball lifetime on the thermal dilepton yield is 
studied (Sec. Ill Cl. Finally, conclusions are drawn and 
an outlook to future studies is given in section m 


II. THE MODEL 

The full theoretical description of the dilepton spec¬ 
tra requires to consider a large number of different pro¬ 
duction processes. At SIS energies, all dileptons stem 
from hadronic sources, in contrast to the situation at 
SPS, RHIC or LHC energies, where a significant contri¬ 
bution is assumed to come from ^-annihilation in the 
quark-gluon plasma [44j . In general one can distinguish 
between two different hadronic contributions, such from 
long-lived particles (especially n and 77 mesons) and those 
from the short-lived light vector mesons (mainly p, but 
also ui and (f>). The former have a life time which is sig¬ 
nificantly larger than the duration of the hot and dense 
stage in a heavy-ion collision. Consequently, almost all 
of the decays of the long-lived mesons into lepton pairs 
will happen in the vacuum, i.e., no modification of the 
spectral shape is expected for those contributions. In 
contrast the latter mesons have a short life time and will 
therefore decay to a significant amount inside the fireball 
and their spectral properties are altered by the medium. 

These differences are also significant for our approach. 
For the 7 r and 77 meson contribution it is neither nec- 
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essary nor adequate to calculate thermal dilepton emis¬ 
sion. The yields of these hadrons are determined directly 
from calculations with the UrQMD transport approach. 
For the p and the u>, however, we calculate the ther¬ 
mal emission from coarse-grained transport simulations 
by application of in-medium spectral functions. Here the 
influence of the space-time evolution of the fireball is im¬ 
manent, as the spectral shape will largely depend on the 
values of temperature T and baryon chemical potential 
Pb- Although for the <\> some medium-modifications of 
the spectral shape are predicted as well, we here skip 
a full thermal calculation for the present investigation. 
On the one hand, <f> production is strongly suppressed 
at the low energies considered here and will therefore 
hardly give any significant contribution to the invariant 
mass spectrum, on the other hand the predicted medium 
effects are rather small. Consequently, the 4> contribu¬ 
tion is directly extracted from the UrQMD calculations, 
as applied for the 7r and rj. At higher invariant dilep¬ 
ton masses, i.e. mainly above 1 GeV/c 2 , also multi-pion 
states in form of broad resonances influence the dilep¬ 
ton production. This contribution is also calculated as 
thermal emission. 

The underlying model for all our considerations is 
the Ultra-relativistic Quantum Molecular Dynamics ap¬ 
proach (UrQMD), which is a non-equilibrium cascade 
model [4511481 . It includes all relevant hadronic reso¬ 
nances up to a mass of 2.2GeV/c 2 . The model gives 
an effective solution to the Boltzmann equation. The 
hadrons are propagated on classical trajectories and can 
interact in form of elastic and inelastic scatterings. Pro¬ 
duction of new particles via resonance formation (e.g., 
7r —{— 7r —> p) or the decay of resonances in form of 
A —> N + 7T. String excitation is possible for hadron- 
hadron collisions with yfs > 3GeV but is negligible in 
the SIS energy regime considered here. 


A. Coarse-Grained Contributions 

For the calculation of thermal emission rates by apply¬ 
ing in-mediunr spectral functions, one needs to extract 
the local thermodynamic properties from the UrQMD 
simulations. To obtain a phase-space distribution that is 
sufficiently smooth in a small volume AU around each 
point in space-time, we simulate a large ensemble of 
events. A grid of space-time cells with Ax = Ay = Az = 
0.7 — 0.8 frn and At = 0.6fm/c is set up and the energy- 
momentum tensor T^ v as well as the net-baryon four-flow 
'j p ‘ are determined in each of these cells as 

J l 1 


The rest frame according to Eckart [TjJj can be found by 
performing a Lorentz boost such that the baryon flow 
vanishes in the cell, i.e. j b = 0. In the Eckart frame we 
extract the energy density e and baryon density pb as 
an input to obtain the temperature T and baryon chem¬ 
ical potential pb by applying an equation of state. At 
SIS energies it is sufficient to use a hadron-resonance gas 
(HG-EoS) [50 • This HG-EoS includes the same degrees 
of freedom as the UrQMD model. 

It is important to bear in mind that the procedure as 
described above assumes local (isotropic) equilibrium in 
the cell. In macroscopic descriptions of heavy-ion col¬ 
lisions kinetic and chemical equilibrium is usually in¬ 
troduced as an ad-hoc assumption. However, we ex¬ 
tract the thermal properties from a transport (i.e. non- 
equilibrium) approach which has effects as viscosity and 
heat conduction. Although the creation of an approxi¬ 
mately equilibrated stage is usually considered to hap¬ 
pen on extremely short time scales after the beginning 
of the collision, it is difficult to prove the creation of 
thermal and chemical equilibrium explicitly. Previous 
studies comparing UrQMD calculations with the results 
from the statistical thermal model showed that it might 
take up to 8-10 fm/c before one can assume the sys¬ 
tem to be in approximate kinetic and chemical equilib¬ 
rium on a global scale ED [52. For practical reasons we 
use the momentum-space anisotropy to characterize to 
which degree the local kinetic equilibrium is constituted 
in the present study. Here, the coarse-grained micro¬ 
scopic transport calculations show a significant deviation 
from equilibrium in form of large pressure differences be¬ 
tween the longitudinal and transverse components of the 
energy momentum tensor at the beginning of the reac¬ 
tion: During the first few fm/c the longitudinal pressure 
is significantly larger than the transverse pressures, which 
is mainly due to the deposition of high longitudinal mo¬ 
menta from the colliding nuclei. It has been suggested 
that in this case one can determine realistic values for the 
energy density by using a generalized equation of state 
[53l [54] where 


£ real — / \ • (2) 

r(x) 

For a system with Boltzmann-type pressure anisotropies 
the relaxation function r( x) takes the form 
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where x = (P||/P_l) 3 ' 4 denotes the pressure anisotropy. 
With this result we can extract meaningful energy densi¬ 
ties from the coarse-grained distributions also in the very 
early stage of the reaction, as has already been shown in 

m- 

With the values of T and pb known for all cells, the 
calculation of the thermal dilepton emission is straight- 
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forward (7[ 153 and takes the form 


d Nu 
d 4 a;d 4 q 


tt 3 M 2 


f B (q-U;T) 


x Imll^(M,q; p B ,T) 


(4) 


with the Bose distribution function / B and the lepton 
phase-space L(M), while M and q denote the invariant 
mass and the momentum of the lepton pair, respectively, 
and a em is the electromagnetic coupling constant. The 
relevant physical quantity is the retarded electromagnetic 
current-current correlator Ilcm^ which contains all the in¬ 
formation on the medium effects of the spectral function. 

Note that equation ([I]) is derived for the case of full 
chemical equilibrium, which requires all chemical poten¬ 
tials of non-conserved charges to be zero. Similar to the 
kinetic anisotropies at the beginning of the reaction, also 
deviations from the chemical equilibrium composition of 
the system appear at the early stages of transport simu¬ 
lations. Especially an overpopulation of pions is observed 
here, which dominates the evolution for a significant pe¬ 
riod of time [56] . The appearance of a finite pion chemical 
potential was explained by the large initial production 
of pions and their long relaxation time [571 HE]- This is of 
particular importance for the present study since a finite 
chemical potential was considered to have a significant 
influence on the thermal dilepton emission rates by its in¬ 
fluence on the 7 T-p interactions j5SJ[6(J- For the present 
calculations, we extract the pion chemical potential in 
Boltzmann approximation for each cell and consider its 
effect on the thermal dilepton emission (i) in form of its 
direct influence on the spectral functions and (ii) as ad¬ 
ditional fugacity factor, which will show up in equation 
© for the chemical non-equilibrium case. 

In the following the different thermal contributions 
which are considered for the present study are discussed. 


1. Thermal p and u> Emission 


butions to the self-energy of the p are taken into account 
here. These are 


(a) the modification of the pion cloud of the p meson by 
particle-hole and A-hole excitations in the medium, 

(b) scattering with mesons (M = 7r, K, K , p) and 

(c) scattering with the most abundant baryon resonances 

(B = N, A1232, Ar* 44 0 ,...). 

The in-medium propagator consequently takes the form 

Dp = M ^ - m2 - - S pM - ' (5) 

The oj spectral function [66] is similarly constructed. 
However, here the situation is more complicated as the 
oj is basically a three-pion state, and the vacuum self en¬ 
ergy is given by decays into pn or 37 t. The full in-medium 
propagator reads 

D u =[M 2 — m 2 u + i m^ (T 37r + T p7r + r w7r _ >7r7r ) 

v v i-i (6) 

^UJTTbi 


It includes the following contributions: 

(a) oj —i pir decays including the corrections from the 
medium-modified p spectral function, 

(b) the direct oj —► 37r decays, 

(c) ojtt —» 7T7T inelastic scattering, 

(d) ojtt — ► bi resonance scattering and 

(e) ojN — > Nn 520 ), IV(i 65 o) resonance scattering at an ef¬ 
fective nucleon density p e ff. 

To take into account the off-equilibrium of the pions an 
additional fugacity factor is introduced in equation Q, as 
already explained above. In one-loop approximation, for 
the thermal p and oj emission one obtains an additional 
overall fugacity factor [7] 


Several approaches for the description of in-medium 
effects on vector mesons exist. However, a full description 
of the different effects that influence the spectral shape 
in a hot and dense medium are highly non-trivial, and 
there are only a few calculations that include both the 
effects of finite temperature and density. E.g. in Ref. [6T] 
the in-medium spectral functions were determined using 
empirical scattering amplitudes. However, this spectral 
function is calculated in low-density approximation for a 
weakly interacting pion-nucleon gas and tested only up 
to densities of 2 p 0 . The application of this approach is 
therefore questionable for the situation at SIS energies 
where very high net baryon densities are reached. 

In the present work a calculation from hadronic many- 
body theory [621465] is applied, which has proven to suc¬ 
cessfully describe the dilepton spectra at SPS and RHIC 
energies [Til GT> . In the medium, three different contri¬ 


■n _ e ^7i„/T 


(7) 


The exponent n hereby takes the value 2 in case of the 
p emission and 3 for the oj contribution |T5] . The final 
yield is then calculated according to 


AN e+e - = A 4 [ &P <FV e +e- r n 

AM J p 0 d 4 a;d 4 p“ 7r ' 


( 8 ) 


Note that the four-volume Aa; 4 = At.AV of the cell is an 
invariant quantity and therefore the same in all reference 
frames. 

As it is not reasonable or possible to calculate thermal 
dilepton emission for all cells (e.g., due to low tempera¬ 
ture) we also include a non-thermal (“freeze-out”) con¬ 
tribution for the p and oo meson. Those dileptons are 
directly extracted from the UrQMD calculations (for de¬ 
tails see section IIB 3). 
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2. Multi-n Contribution 

Although it is known from previous theoretical in¬ 
vestigations at higher collision energies that the multi- 
pion contribution plays a significant role only for masses 
greater than 1 GeV/c 2 , which is mostly beyond the kine¬ 
matic limit of dilepton production at SIS energies, we in¬ 
clude this contribution for completeness. Here the same 
description as developed in mnu is applied which uses 
chiral reduction techniques for the case of the chiral limit 
EZ]. This leads to a chiral mixing of the isovector part 
of the vector and axial-vector correlators. The isovector- 
vector current correlation function takes the form 

n v (p) =(i - e)zi f 

£ ( 9 ) 

+ 2 (4+ 4) 

where e denotes the mixing parameter which depends on 
temperature, critical temperature and the pion chemical 
potential. The result is an admixture of three-pion and 
five-pion axial-vector pieces on the vector four-pion part. 
Once again, as in case of the p and w spectral functions, 
the effect of a finite pion chemical potential also enters 
in form of a fugacity factor z" with n = 3,4, 5 denoting 
the pion multiplicity of the corresponding state. Note 
that the two-pion and the three-pion piece corresponding 
to the decay ai — > p + it are already included in the p 
spectral function and are therefore not considered for the 
multi-pion contribution. 


B. Transport Contributions 

1. 7r and rj Dalitz Decay 

The 7 r° and the rj pseudo-scalar mesons can both decay 
into a lepton pair via the Dalitz decays 

7T° —>■ 7 l + l~, (10) 

Following the scheme presented in we treat this de¬ 
cay as a two-step process: First the decay of the pseudo¬ 
scalar meson into a photon 7 and a virtual photon 7 * 
and the following electromagnetic conversion of the 7 * 
into a lepton pair [BS 1 . The width of the meson decaying 
into a photon and a virtual photon can be related to the 
radiative decay width. The full expression then takes the 
form 


agreement with the values as predicted by the Vector 
Dominance Model [ 68 l : 


f> 77 »(m 2 ) 

f„ 77 .(m 2 ) 


1 + b^o IV/ 2 , 



( 12 ) 


with 0 = 5.5 GeV -2 and A v = 0.72 GeV. Note that for 
this contribution, we only consider those particles from 
the final freeze-out of the calculation, neglecting all n 
and 77 which are reabsorbed during the evolution of the 
collision. 


2. (j> Direct Decay 

The width for the direct decays of the vector mesons 
can be determined via [551 

!WM) = Tv ^™ v) v ±L{M). (13) 

For the branching ratio at the meson pole mass 
my we use the value from the particle data group, 

r^ ee (m v )/r tot = 2.95 • 10- 4 [7DJ. 

An important difference compared to the procedure 
for the long-lived pseudo-scalar mesons is the assump¬ 
tion that the (f> mesons continuously emit dileptons over 
their lifetime ED- That means, we track the time of 
production and the decay (or absorption) of the 4> and 
integrate over the corresponding lifetime in the particle’s 
rest frame. 

d N lt A N u ^ ft dt T v ^ u (M) 

AM AM ^ ^ J t iv AM ' ^ ’ 

The factor yy is here introduced to account for the rel¬ 
ativistic time dilation in the computational frame com¬ 
pared to the vector mesons rest frame. 

This “shining approach” is identical with the assump¬ 
tion of simply having one additional dilepton from each 
<j> (weighted with the according branching ratio) if there 
is no absorption. But inside a dense medium, there is 
a significant chance that the particle will not decay but 
suffer an inelastic collision. Therefore the probability for 
the decay into a dilepton is reduced, and the shining as 
described above accounts for this effect m- 


3. Non-thermal p and to 


dr 


M-tjl+l- 

d M 2 


= 2 r 


M-¥ 27 


1 - 


M 2 


m 


M 


( 11 ) 


X |Fm 77 .(M 2 )P 


3ttM 2 


L(M). 


With the lepton phase space L(M) and the form factor 
Fp 77 . (M 2 ). Here we use the form factors as obtained 
from fits to experimental data, which are in very good 


There are two situations, where a thermal calculation 
of the dilepton emission from spectral functions becomes 
difficult or even unreasonable. Firstly, for cells with no 
baryon content, where the Eckart definition of the rest- 
frame does not apply (those cells are found in the late 
stage of the reaction in very peripheral cells) and sec¬ 
ondly for cells where the temperature is found to be be¬ 
low 50 MeV. In the latter case, when going to such usually 
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Figure 1. (Color online) Time evolution of the energy and baryon densities e and p in units of the according ground-state 
densities (a) and for temperature T as well as the baryon and pion chemical potential, /xb respectively p n (b) for the central 
cell (i.e., at x = y = z = 0). The results as obtained via coarse-graining of microscopic transport calculations are shown for 
Ar+KCl collisions at _Ei a b = 1.76 AGeV and for Au+Au reactions at 1.23 AGeV. 


low-density cells, the determination of T and pb becomes 
less accurate and one will necessarily come to the point 
where the assumption of a thermalized system in the cell 
becomes unreliable. Consequently, we do not assume any 
thermal dilepton emission here. This procedure is in line 
with the findings of thermal-model studies d, where it 
was shown that the freeze-out temperature in heavy-ion 
collisions at E\ a \> = 1 — 2 AGeV is around 50 MeV. This 
also indicates that it is neither necessary nor suggestive 
to assume thermalization of the system at lower temper¬ 
atures. 

However, one has to consider that dilepton emission 
from p and <x> mesons is of course also possible in the 
cells for which one of the conditions mentioned above per¬ 
tains. As the macroscopic picture is questionable here, 
we apply a similar procedure as for the <f> (described in 
section IIB 2) to extract the clilepton emission from the 


microscopic simulation. The width for the direct decays 
of the vector mesons can likewise be determined via equa¬ 


tion (13). For the branching ratio at the meson-pole 


mass we use the values from the particle data group, i.e. 


p—>ee 


/r tot = 4.72 • 10- 5 and r^ ee /T tot = 7.28 • 10 


-5 


HOI- However, the u> can not only decay into a lepton 
pair directly, but also via Dalitz conversion into a pseudo¬ 
scalar meson and a dilepton. As in the case of the n and 


the r] equation (11) applies here. Only the form factor is 
different and takes the form 


\F U (M*)\* = 


A 2 (A 2 

CO \ Cl 


■77 


(A 2 - M 2 ) + A 2 y 2 


(15) 


with the parameters A w = 0.65 GeV and 7 ^ = 0.04 GeV 

mm- 

Similarly to the procedure for the </> meson, a continu¬ 
ous emission of dileptons from the p and u> is assumed in 
these special cases. However, as we consider space-time 



Figure 2. (Color online) Time evolution of the thermal vol¬ 
ume V(p M ' as seen from the center-of-momentum frame of 
the collision. This is equal to the thermal four-voulume 
V4 for each timestep divided by the length of At, which is 
0.6fm/c here. The results are shown for Ar+KCl collisions 
at Ai a b = 1.76 AGeV and for Au+Au reactions at 1.23 AGeV 
as obtained via coarse-graining of microscopic transport cal¬ 
culations. The results are plotted for different temperatures. 


cells with a definite length of the time-steps At, the dilep¬ 
ton rate is multiplied with this time instead of the actual 
particle’s lifetime within the cascade simulation. This 
is done to guarantee consistency, avoid double counting 
and strictly distinguish between cells with thermal and 
non-thermal emission. Consequently, here equation (14) 


















7 




Figure 3. (Color online) Invariant mass spectra of the dielectron yield for Ar+KCl collisions at E\ a b = 1.76 AGeV (a) and for 
Au+Au at Slab = 1.23 AGeV (b). The results are normalized to the average total number of 7r° per event and shown within 
the HADES acceptance. The results for Ar+KCl are compared to the experimental data from the HADES Collaboration [23]. 


takes the form 

d Nu _ ANu _ AtT v ^ n {M) 

AM AM 1 AM 

i— 1 j —1 lv 

and is applied for each time step At. 


III. RESULTS 

For the results presented here we used calculations 
with an ensemble of 1000 UrQMD events. However, sev¬ 
eral runs using different UrQMD events as input had to 
be performed to obtain enough statistics especially for 
the non-thermal p and ui contributions. Note that in 
case of the experimental Ar+KCl reaction we simulated 
the collision of two calcium ions instead, as this makes 
the calculation easier for symmetry reasons. Effectively 
it is the same as the Ar+K or Ar+Cl reactions that were 
measured in the experiment and the size of the system 
remains identical. To simulate the correct impact pa¬ 
rameter distribution, we made a Woods-Saxon type fit 
to the HADES trigger conditions for Ar+KCl [75 and 
Au+Au m- In both cases this approximately corre¬ 
sponds to a selection of the 0-40% most central colli¬ 
sions. The number of n° per event, which will be impor¬ 
tant for the normalization of the dilepton spectra, are 
found to agree well with the HADES measurement for 
Ar+KCl reactions. Here the HADES collaboration mea¬ 
sured N°o P = 3.5 where we find iV sl 0 m « 3.9, i.e. the 
deviation is only 12%. For the larger Au+Au system 
a number iV® 1 ™ ss 8.0 results from the events generated 
with the UrQMD model. Note that for reasons of self- 
consistency we normalize the dilepton spectra with the 
UrQMD 7r° yield, not the experimental one. 


The final dilepton results were filtered with the 
HADES acceptance filter EH, and momentum cuts were 
applied to compare the simulations with the experimen¬ 
tal results. As only very preliminary results and no filters 
are available for the Au+Au reactions at 1.23 AGeV, we 
used the same filter as for p+p and p+n reactions at 
1.25 AGeV which should be quite close to the final ac¬ 
ceptance [761. 

In case of the DLS Ca+Ca spectrum, version 4.1 of the 
DLS acceptance filter m is used. Furthermore, an RMS 
smearing of 10% is applied to account for the detector 
resolution. For this reaction we used a minimum-bias 
simulation of Ca+Ca events, because impact-parameter 
distributions are not available for DLS. Here the final 
invariant-mass spectrum is normalized to the total cross- 
section of a Ca+Ca reaction. 


A. Reaction Dynamics 

The main difference between the two heavy-ion reac¬ 
tions considered here (as measured by the HADES Col¬ 
laboration) is the size of the colliding nuclei. Therefore, 
it is interesting to first have a look at the evolution of 
the reaction for both systems. In Fig. [l] the evolution of 
the baryon and energy density (a) as well as the evolu¬ 
tion of temperature and the chemical potentials (b) are 
shown. The maximum density values in the central cell 
of the grid, i.e., in the center of the collision, reach sim¬ 
ilar values up to roughly 3-6 times ground-state baryon 
density po and energy density £+ for both reactions. In 
case of the larger system (Au+Au) a plateau develops for 
a duration of more than 10 fm/c, while for Ar+KCl (re¬ 
spectively Ca+Ca in our simulations) the densities drop 
off rather quickly after reaching the maximum. Note that 
the values for the energy density e shown in Fig.[l](a) are 
















Figure 4. (Color online) Comparison of invariant mass spectra with the full spectral function and for the case of no baryonic 
effects (i.e. for p e ff = 0). As in Fig. [ 3 } the dielectron yields for Ar+KCl collisions at i?i a b = 1.76 AGeV (a) and for Au+Au 
at £i a b = 1.23 AGeV (b) are shown within HADES acceptance and normalized to the average number of produced n°. Note 
that the UrQMD contributions are included in the sum, but the different single yields are not shown explicitly for reasons of 
lucidity. 


corrected for the pressure anisotropy at the beginning of 
the collision. 

With regard to the evolution of temperature and bary- 
ochemical potential in Fig. l](b) , we find that again sim¬ 
ilar top values are obtained. In both reactions, Ar+KCl 
and Au+Au, we get peak values of T ss 100 MeV and 
Pb ~ 900 MeV. (Note that the quark chemical potential 
p q = |/iB is shown instead of the baryon chemical po¬ 
tential.) Especially for the Au+Au reaction the baryon 
chemical potential shows a more prominent plateau than 
the baryon density. However, note that T and pb depend 
on e and pb non-linearly via the EoS. For the collisions 
at SIS energies studied here, pb rapidly rises to values 
very close to the nucleon mass - but once it reaches such 
a high level it exhibits a much less significant rise in spite 
of a further increase of pb- This is a consequence of the 
Fermi statistics which takes effect in this case. Further¬ 
more, it is clear that the values of pb are much higher 
here and show a different evolution than in our recent 
study for top SPS energy [TB]. In case of Au+Au col¬ 
lisions the central cell stays for approximately 20fm/c 
in a stage with extremely high baryochemical potential, 
so that any baryon-driven effect on the dilepton spectra 
should be clearly visible for this system. Similar findings 
are also true for Ar+KCl reactions, but with significantly 
shorter lifetime. The maximum temperature is approxi¬ 
mately 10 MeV higher in the latter reaction due to the 
slightly higher collision energy. 

The pion chemical potential rises up to values of 
around 100 MeV and then equally quickly drops to val¬ 
ues around 20-50 MeV for the rest of the reaction. The 
peak at the beginning can be explained with the non¬ 
equilibrium nature of the cascade, where a large number 
of pions is produced rapidly at the very beginning before 


the system has time to equilibrate. 

The corresponding evolution of the thermal four- 
volume at each time step, divided by its duration At, 
as obtained for both systems from the coarse-grained mi¬ 
croscopic transport calculations is shown in Fig. [2] This 
quantity is the spatial volume for the different tempera¬ 
ture ranges as seen from the center-of-momentum frame 
of the collision. One observes that the hotter cells cre¬ 
ated during the evolution of the system (for tempera¬ 
tures above 80 respectively 90 MeV) are present only for 
comparatively short time spans, and their number is sig¬ 
nificantly smaller than the volume of cells with lower 
temperatures. While we only show the evolution of the 
first 30fm/c in Fig. [ 2 ] for the sake of clarity, we men¬ 
tion that even at the end of the simulation after 70fm/c 
one still finds a few cells at temperatures above 50 MeV. 
Note, however, that the dilepton emission from these late 
stage cells is marginal and insignificant with respect to 
the overall time-integrated dilepton spectra. 


B. Dilepton Spectra 

1. HADES Results 

The resulting dielectron invariant-mass spectra for the 
two heavy-ion reactions measured by the HADES Col¬ 
laboration are presented in Fig. [3j The left figure (a) 
shows the Ar+KCl results at 1.76 AGeV compared to 
the experimental data and the right figure (b) shows our 
prediction for the larger Au+Au system at 1.23 AGeV 
(no dilepton data are published for the Au+Au measure¬ 
ment yet). The coarse-graining results for the Ar+KCl 
case show good agreement with the experimental data, 
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Figure 5. (Color online) Transverse-mass spectra of the dielectron yield in different mass bins (a)-(e) and the dielectron rapidity 
spectrum of pairs with invariant mass M ee < 0.13 GeV/c 2 (f) for Ar+KCl collisions at i?i a b = 1.76 AGeV. The results are 
compared to the experimental data from the HADES Collaboration [25]. In contrast to Figures [3] and [4] the theoretical results 
here are not shown within the HADES acceptance, as these data are already fully corrected for acceptance and efficiency. 


especially the mass range between 0.2 and 0.6 GeV/c 2 
can be well described within this approach. This is dif¬ 
ferent from previous simulations which have indicated an 
excess of the experimental outcome above the cocktail 
respectively the transport calculations 23 79, 201. The 
dominant contribution stems from a broadened p. It is 
further noteworthy that also the ui shows a non-negligible 
broadening in our model results. A slight overestimation 
of the experimental dilepton yield shows up for the low- 


mass region (below 150 MeV/c 2 ) which is dominated by n 
Dalitz decays, an effect that also has been found in other 
transport calculations E3- However, it remains unclear 
where this pion excess stems from. As the dielectron yield 
is normalized to the total ir° number, theory and exper¬ 
iment should agree in the pion-dominated mass-region. 
Therefore the deviation may be due to a phase-space ef¬ 
fect, connected to the geometrical detector acceptance. 
Another slight excess of the model results is manifest 
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around the pole mass of the p. A significant part of the 
dileptons here stems from non-thermal p mesons, which 
directly come from the UrQMD calculations using the 
shining method. As for UrQMD the p-production cross 
sections in the threshold region are known to be slightly 
too large |S 0 j, this overestimation is probably due to the 
non-thermal contribution. 

The medium effects become even stronger in the larger 
Au+Au system, as can be seen in the right part (b) of 
Fig-! The yield from the thermal p is higher at low 
masses compared to the Ar+KCl reaction and also the 
thermal ui is slightly stronger here. This is in line with 
the findings of the previous section MB where it was 
pointed out that the hot and dense stage lives substan¬ 
tially longer in the Au+Au reaction. Consequently, the p 
contribution is larger for the Au+Au system as compared 
to the Ar+KCl reaction, especially at the low masses 
where the in-medium broadening comes into account. 

The dramatic effect of the presence of baryonic matter 
on the dilepton spectra can be seen in Fig. |4j Here the 
thermal contributions to the e + e _ invariant mass spec¬ 
trum of the p and the w for the case that no baryons and 
anti-baryons are present (i.e. the effective baryon den¬ 
sity Peff is put to zero) are compared to the results with 
the full medium effects. The baryon-driven effects are 
significant, with the broadening around the pole masses 
of the p and ui meson and the strong increase in the low- 
mass dilepton yield, which is again more distinct for the 
larger Au+Au system as compared to Ar+KCl. Addi¬ 
tionally the total sum for the two cases with and without 
baryonic effects is shown, including also all the UrQMD 
contributions (details are omitted for clarity). Here we 
find stronger differences between Ar+KCl and Au+Au, 
mainly due to the relatively smaller contribution from 
the 77 for Au+Au reactions, compared to the stronger 
broadening of the p. 

Looking not only at the invariant-mass spectra but also 
at the transverse-mass distributions in different mass bins 
in Fig.[5](a)-(e), one finds again a good agreement of our 
model calculations with the HADES data for Ar+KCl. 
While for the lowest mass bin (M ee < 0.13GeV/c 2 ) the 
pion contribution dominates, the thermal p is the most 
significant contribution to the dileptons in all other bins. 
For the highest mass bin ( M ee > 0.65GeV/c 2 ), which 
includes the pole mass region of the p, one can observe a 
slight overestimation of the total yield, similar to the one 
observed in the invariant mass spectrum in this mass re¬ 
gion (see Fig.[3](a)). Once again, we argue that this is due 
to the high p-production cross section in the threshold re¬ 
gion and therefore stems from the non-thermal transport 
p. As one can see, the thermal p alone would describe 
the data very well. 

Note that in some kinematic regions statistical fluc¬ 
tuations are seen in the transverse mass spectra. How¬ 
ever, naturally this mainly affects subleading contribu¬ 
tions and the very high transverse masses. As the pro¬ 
duction of certain particles is highly suppressed at SIS 
energies (e.g., in case of the <j) meson) it would need ex- 



Figure 6. (Color online) Invariant-mass spectrum of the di¬ 
electron yield for Ca+Ca collisions at I?i a b = 1.04 AGeV 
within the experimental acceptance. The result is compared 
to the data from the DLS Collaboration |19| . 

cessive computing resources to remove these statistical 
fluctuations completely. Nevertheless, the effect on the 
total yield is usually rather small - especially compared 
to the uncertainty of the experimental measurement. In 
general, an estimate of the global (statistical and system¬ 
atic) error of our dilepton calculations is extremely diffi¬ 
cult due to the many different parameters (cross-sections, 
branching ratios, spectral functions, etc.) which enter the 
calculation and because for other sources (e.g., the filter¬ 
ing) the error can hardly be quantified. At least for the 
long-lived contributions (mainly 7 r and p) a comparison 
between different transport models indicates that their 
contribution to the dilepton spectra is quite well deter¬ 
mined and does not differ much between the models w 
[22]. For the thermal contributions (p and w) the two 
main error sources are the uncertainties of the spectral 
function and of the description of the reaction dynam¬ 
ics, i.e., the time-evolution of temperature and chemical 
potential. 

For the lowest mass bin, the HADES collaboration has 
also measured the rapidity distribution of dielectrons. 
The results from the coarse-grained simulations as well as 
the data points are shown in Fig. |5](f). The shape of the 
spectrum is reproduced well, but with some 20 % excess 
above the data points which was already visible in the 
invariant mass spectrum. Note that the rapidity values 
are for the laboratory frame, i.e. for £i a b = 1.76 AGeV 
mid-rapidity corresponds to a value of yo = 0 . 86 . 


2. DLS Results 

Although the DLS measurement was done with a 
smaller acceptance and lower resolution than in the more 
recent analyses by the HADES Collaboration, it is nev- 
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Figure 7. (Color online) (a) Ratio of the thermal four-volume V4 for different temperatures to the mass number A of the 
colliding nuclei. The results are normalized to the ratio obtained with 12 C + 12 C collisions, (b) Time duration over which the 
central cell of the coarse-graining grid (for x = y = z = 0 ) emits thermal dileptons, i.e., for this period of time the central cell 
has a temperature which is above 50 MeV (blue squares) and the same scaled by A 1 / 3 (green triangles). Both plots (a) and 
(b) show the results in dependence on mass number A for central collisions and a collision energy = 1.76 AGeV. 


ertheless interesting to compare our results also with 
these heavy-ion data. The system Ca+Ca which was 
measured by DLS is comparable to the Ar+KCl re¬ 
action, as mentioned above. However, the DLS Col¬ 
laboration performed the measurement at a lower en¬ 
ergy of £?i a b = 1.04 AGeV. In Fig. [6] the coarse-grained 
UrQMD results for the according invariant-mass spec¬ 
trum is shown within the DLS acceptance, together with 
the experimental data points m- In general the spec¬ 
trum does not differ strongly from the simulations for 
HADES, but the peaks - especially in case of the u> - 
are more smeared out here, which is due to the low mass 
resolution of only cxm/AI = 10 % of the detector, and 
the acceptance at low masses is suppressed compared to 
the HADES measurements. (For comparison, the mass 
resolution of the HADES experiment is roughly 2% in 
the region of the p and ui pole masses [30].) While the 
overall description of the experimental spectrum with the 
present model is relatively good, a slight excess of the 
data above our model curve is present in the mass range 
from 0.2 up to 0.6 GeV/c 2 . This is in contrast to the 
findings for the HADES case, where we could describe 
this mass region quite accurately. 

What might be the reasons for this deviation? It is 
possible that at this very low energy other processes 
presently not considered in our model might become 
more dominant here, e.g., an explicit bremsstrahlung 
contribution (note, however, that some bremsstrahlung 
effects are already considered within the in-medium spec¬ 
tral functions). As it has been shown that the importance 
of bremsstrahlung contributions increases with decreas¬ 
ing collision energy, these effects might be more signif¬ 
icant than in Ar+KCl reactions at the slightly higher 
energies used for the HADES experiment. Nevertheless, 
there are still uncertainties within the different model 


calculations for the bremsstrahlung contribution, differ¬ 
ing by a factor 2-4 [HU 1E2- A main problem here is 
the correct determination of the overall effect of the dif¬ 
ferent interfering processes, which is highly non-trivial. 
Another issue is that the lower collision energy also re¬ 
sults in slightly lower temperatures peaking around 80 
MeV, which is indicating less thermalization of the col¬ 
liding system. Besides, one has to take into account that 
the DLS experiment with its two-arm set-up had a more 
limited acceptance as compared to the HADES detector. 
This makes it difficult to draw clear quantitative conclu¬ 
sions from the comparison to the data. Additionally, the 
lack of a measured impact-parameter distribution ham¬ 
pers precise calculations within the coarse-graining ap¬ 
proach, as the medium effects can be quite sensitive to 
the centrality of the collision. 

However, it was shown (at least for C+C reactions) 
that the HADES results agree very well with those mea¬ 
sured by the DLS collaboration if the HADES data are fil¬ 
tered with the DLS acceptance [21] . So these data should 
provide a good additional check for theoretical models in 
spite of their lower accuracy. Overall, the agreement of 
our results with the data is quite good, especially consid¬ 
ering that previous pure transport calculations with the 
UrQMD model (without bremsstrahlung) clearly under¬ 
estimated the Ca+Ca data from DLS by a factor of 5-10 
in the mass range from 0.2 to 0.6 GeV/c 2 [72] . 

C. System size and lifetime effects on thermal 
dilepton emission 

The results of the previous section indicate that the 
size of the colliding nuclei (and in consequence also the 
size and duration of the hot and dense fireball created 
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Figure 8. (Color online) (a) Ratio of the thermal (red squares) and non-thermal dilepton yield (blue triangles) in the invariant 
mass range from 0.2 to 0.4GeV/c 2 to the mass number A of the colliding nuclei; the same ratio is plotted for the total number 
of 7 r° (green triangles). The results are normalized to the ratio obtained with 12 C + 12 C collisions, (b) Ratio of the thermal 
dilepton yield in the invariant mass range from 0.2 to 0.4GeV/c 2 to A 4 ^ 3 (red squares), to the thermal four-volume (blue 
triangles), to the product of A and the time f thermal in which the central cell (for x = y = z = 0) emits thermal dileptons, i.e, 
is at T > 50MeV (green triangles), and to the number of 7r° produced in the reaction scaled with an exponent 4/3 (orange 
circles). All ratios are normalized to the result obtained for 12 C + 12 C. Both plots (a) and (b) show the results for central 
collisions and a collision energy = 1.76 AGeV. 



Figure 9. (Color online) Ratio of the thermal to the non- 
thermal dilepton yield in the mass range 0.2 GeV/c 2 < M ee < 
0.4GeV/c 2 for A+A collisions at E\^h = 1.76 AGeV. This 
ratio depicts the “excess” of the thermal dilepton yield with 
regard to the cocktail contributions. The results are shown in 
dependence on the mass number A of the colliding nuclei. 


thereby) largely influence the thermal dilepton yield. 
This encourages a systematic study of the system size de¬ 
pendence for thermal and non-thermal contributions and 
also raises the question, up to which point the assumption 
of a thermalized system seems reasonable. For this pur¬ 
pose we compare different systems from C+C to Au+Au 
in central collisions at an energy of E\ & \, = 1.76 AGeV. 


Fig. [T](a) shows the ratio of the thermal four-volume V 4 
for different temperatures to the mass number A of the 
colliding nuclei. The results are normalized to the ratio 
obtained with 12 C+ 12 C collisions. It is noteworthy here 
that the four-volume of the hottest cells (T > 90MeV) 
shows a much stronger increase than the overall thermal 
volume (i.e., for T > 50MeV). However, only the rela¬ 
tive increase is shown in this plot. In direct comparison 
the number of cells with highest temperatures is very 
small compared to the total volume (approximately 1 /30 
for C+C). It is furthermore striking that the increase of 
the thermal four-volume is not proportional to A, but it 
shows a stronger increase for larger systems. This is not 
surprising, as A is only a measure for the volume of the 
colliding nuclei. In studies with the statistical model it 
was shown that at SIS energies the thermal freeze-out 
volume is closely related to the initial overlap volume of 
the system created in A — A collisions [73] [83]. However, 
the thermal four-volume is also determined by the life¬ 
time of the fireball. As it is difficult to determine some 
kind of an average lifetime within the coarse-graining ap¬ 
proach, we concentrate on the central cell of the grid and 
assume that the time for which it remains at an tempera¬ 
ture greater than 50 MeV should be to first order a good 
approximation of the overall thermal lifetime. Fig. [7](b) 
shows the mentioned time duration over which the cen¬ 
tral cell emits thermal dileptons in dependence on A. 
Note that this duration approximately scales with A 1 / 3 , 
i.e. with the diameter of the nuclei. Obviously the main 
influence on the lifetime at those low energies seems to 
be the time the two nuclei need to traverse each other. 

The question now is which influence the evolution of 
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the four-volume has on the production of lepton pairs. 
In Fig. |8](a) the ratio of thermal and non-thermal dilep¬ 
ton yield as well as the total 7 r° yield in the invariant 
mass range from 0.2 to 0.4 GeV/c 2 in relation to the mass 
number A of the colliding nuclei is shown. The results 
are - again - normalized to the ratio obtained with C+C 
collisions. While the ratio of the 7r° number as well as 
the non-thermal dileptons to the mass number remains 
roughly 1 for all system sizes (i.e. the non-thermal dilep¬ 
ton and 7T° yield increases linearly with A ), the thermal 
yield shows a significantly stronger rise. This finding is 
in line with the larger thermal excess found in our calcu¬ 
lations in Au+Au compared to Ar+KCl reactions. 

If one now compares the thermal yield with the cor¬ 
responding total thermal four-volume (which is all cells 
with T > 50MeV), as shown in Fig. [8](b) , the finding 
indicates that the ratio between both shows only a very 
slight increase and remains almost constant, independent 
of the system size. The same is found if one calculates 
the relation between the thermal yield and the product 
of A with the thermal lifetime of the central cell (com¬ 
pare Fig.[7](b)) or directly looks at the ratio between 
the thermal dilepton yield and mass number A respec¬ 
tively the number of produced neutral pions N n 0 scaled 
by an exponent a = 4/3. All those quantities give an 
approximate measure of the thermal four-volume at the 
low energies considered here. At higher collision ener¬ 
gies as, e.g., obtained at the CERN SPS or at RHIC 
(where the whole fireball is more pion-dominated rather 
than baryon-dominated) one would still expect an in¬ 
crease with Af“ 0 , but not any longer with A a . However, 
note that in these cases the nuclei will traverse each other 
much faster and one will also find a significant transverse 
expansion of the fireball, so that the diameter of the nu¬ 
clei can no longer be considered as a rough measure of the 
lifetime; consequently the parameter a might be different 
here. Nevertheless, in a different model calculation [33] it 
was found that the thermal dilepton yield at top RHIC 
energy scales with the number of charged particles as 
N " h , where a is found to take a value of approximately 
1.45. This result is not so far from the value obtained 
within our simple picture - at a completely different en¬ 
ergy regime. 

In consequence, we learn from these results that the 
non-thermal dilepton contributions (from long-lived me¬ 
sons as the 7T° and the 77 , respectively from the freeze-out 
p and uj) increase with A and therefore directly with the 
volume of the colliding system. This is due to the fact 
that these contributions reflect the final hadronic com¬ 
position, after the whole reaction dynamics has ended. 
Therefore the time-evolution of the system is irrelevant 
here. In contrast, the thermal dilepton emission does not 
only increase with the volume, but also with the time in 
which the colliding system remains in a hot and dense 
stage. This is obvious, as the thermally emitted dilep¬ 
tons will escape the fireball unscathed, while all hadronic 
particles undergo processes as rescattering and reabsorp¬ 
tion. In consequence, the thermal yield increases roughly 


proportional to V ■ At ~ A • A 1 / 3 = A 4 / 3 . Exactly due 
to the different mechanisms contributing to thermal and 
non-thermal dilepton yields, one observes a significant 
increase of their ratio when going to larger system sizes. 
This ratio is shown (for the invariant mass range form 
0.2 to 0.4GeV/c 2 ) in Fig. [ 9 ] with is value increasing from 
slightly over two for C+C up to 14 for Au+Au. 

The present results qualitatively agree with a recent 
study on the same issue performed within different mi¬ 
croscopic transport models, where the total dilepton yield 
was found not to scale with A or the number of neutral 
pions, but showing a stronger increase due to the com¬ 
plicated dynamics of the reaction [25]. There, from a 
microscopic point of view, it was also argued that the 
time evolution of the reaction is the main reason for 
the increasing dilepton yield in larger systems. If the 
dense phase with many binary scatterings lasts longer, 
according to this picture, a larger bremsstrahlung contri¬ 
bution (which is proportional to the number of collisions) 
and the repeated regeneration of A resonances raises the 
dilepton yield. However, the contribution of other bary- 
onic resonances was not explicitly considered there. 


IV. CONCLUSIONS & OUTLOOK 

We have presented results on dilepton production in 
Ar+KCl and Au+Au collisions at GSI SIS 18 and in 
Ca+Ca collisions at BEVALAC energies. The results 
are obtained using a coarse-grained microscopic trans¬ 
port approach and employing state-of-the-art spectral 
functions. With this approach the experimental dilepton 
spectra in heavy-ion collisions at E\^ = l-2AGeV can 
be successfully described. The model represents a third 
way to explore the dynamics of heavy-ion reactions. In 
contrast to hydrodynamic/thermal fireball calculations 
or microscopic transport simulations it allows for a con¬ 
sistent treatment of both, high-energy and low-energy 
collisions. 

Our results show that the dominant in-medium effect, 
which is also visible in the dilepton spectra, stems from 
a strong broadening of the spectral shape of the p me¬ 
son due to the high baryon density in the fireball. This 
causes a melting of the p peak and results in an enormous 
increase of the dilepton production below the pole mass. 
The reason for this is mainly the interaction with the 
baryonic resonances, especially the A and -/V/ 520 , which 
are included in the spectral function and give significant 
additional strength to the p contribution at these low 
masses. The effect is much stronger at SIS energies than 
for RHIC or SPS. Here, the baryon-chemical potential 
remains very high for a significant time. This is visible 
in the dilepton spectra which always represent a four- 
volume integral over the whole space-time evolution. In 
the present study we also find a significant broadening of 
the uj meson at SIS energies. 

Furthermore, the influence of the size of the colliding 
system on the thermal dilepton yield was studied. While 
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we find the non-thermal contribution in the mass range 
from 0.2 to 0.4GeV/c 2 (which stems mainly from the 
long-lived 77 -meson but also from the “freeze-out” p and 
w) to scale linearly with the mass number A, a stronger 
increase of the thermal p and u> yield is observed. Their 
contribution to the dilepton yield scales with yR thermal , 
i.e., the system’s volume multiplied with the lifetime of 
the thermal stage. Since t thermal increases with the di¬ 
ameter of the nuclei which is proportional to A 1 / 3 , we ar¬ 
gue that the thermal dilepton emission should scale with 
A 4 5 6 7 8 9 10 11 / 3 . As the number of neutral pions directly scales with 
A, we find the thermal dilepton yield also increasing with 
N^ 0 3 . For future studies it might be quite instructive to 
check whether this proportionality still holds for higher 
collision energies. 

It is interesting to compare our findings with the results 
of pure transport calculations. There, in recent calcu¬ 
lations the dilepton excess in the invariant-mass spectra 
above the cocktail were mainly explained by two different 
effects: (i) Bremsstrahlung contributions and (ii) Dalitz 
decays of baryonic resonances j26j [29]. It is important 
to understand that both effects are also included in the 
spectral functions applied here. These processes, as im¬ 
plemented in the transport models, correspond to cuts of 
the in-medium self-energy diagrams of the p meson. For 
example, a contribution to the self-energy from a A-hole 
excitation would be represented in the transport model 
by the process A —> pN —> 7 *N, assuming strict vector 
meson dominance. However, due to the completely differ¬ 
ent character of the approaches it is difficult to compare 
the single contributions quantitatively. Also note that 
the self-energies in hadronic many-body quantum-field 
theory approach correspond to a coherent superposition 
of the various scattering processes, while in transport ap¬ 
proaches naturally the various processes are summed in¬ 
coherently. A future detailed analysis of the self-energy 
contributions and their relative strengths might be fruit¬ 
ful to better understand how far the both approaches, 
i.e., the transport (kinetic) and the equilibrium thermal 
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